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We derive nonstatic spherically symmetric solutions of a null fluid, in five dimension (5D), to 
Einstein- Yang-Mills (EYM) equations with the coupling of Gauss-Bonnet (GB) combination of 
quadratic curvature terms, namely, 5D-EYMGB radiating black hole solution. It is shown that, 
in the limit, we can recover known radiating black hole solutions. The spherically symmetric known 
5D static black hole solutions are also retrieved. The effect of the GB term and Yang-Mills (YM) 
gauge charge on the structure and location of horizons, of the 5D radiating black hole, is also 
discussed. 
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I. INTRODUCTION 



Recent years have witnessed a renewed interest to 
study black hole solutions in string-generated gravity 
models which mainly is accomplished by studying so- 
lutions of the Einstein theory supplemented by Gauss- 
Bonnet (GB) term [ll, String theory also predicts 
quantum corrections to classical gravity theory and the 
GB term is the only one leading to second order differ- 
ential equations in the metric. On the other hand the 
black hole solutions in gravity coupled to fields of differ- 
ent types have always drew in much attention, in par- 
ticular, a great interest in solutions to Einstein- Yang- 
Mills (EYM) systems Wu and Yang [HI ob- 
tained static symmetric solution of Yang-Mills equation 
for the isospin gauge group S0{3). The remarkable fea- 
ture of this Wu-Yang ansatz is that the field has no con- 
tribution from gradient and instead has pure YM non- 
Abelian component. A curved-space generalization of 
the Wu-Yang solutions 11] for the gauge group 50(3) 
is shown to be a special case of Yasskin's [31 solutions. 
It is known that non-Abelian gauge theory coupled to 
gravitation, i.e., EYM results to precisely the geometry 
of Reissner-Nordstrom with the charge that determines 
the geometry is gauge charge Indeed, Yasskin (sj 
gave an explicit theorem so that from each solution of the 
Einstein-Maxwell equations one can get solutions of EYM 
equations. One would like to study how these features 
get modified in higher-dimensional (HD) spacetimes and 
whether this theorem holds in HD spacetimes. Recent de- 
velopments in string theory indicate that gravity may be 
truly HD theory, becoming effectively four-dimensional 
(4D) at lower energies. Since non-Abelian gauge fields 
also feature in the low energy effective action of string 
theory, it is interesting to study the properties of the 
corresponding EYM inpresence of GB terms. Mazhari- 
mousavi and Halilsoy [a, 01 have found static spherically 
symmetric HD black hole solutions to coupled set of equa- 
tions of the EYMGB, for SO{N~l) gauge group, systems 
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which are based on the Wu and Yang [Tlj ansatz. The 
corresponding static topological black holes have been 
found independently by others (8l-[l0|. 

It would be interesting to further consider nonstatic 
generalization of Mazharimousavi and Halilsoy solutions 
[a, 0| • It is the purpose of this Letter to obtain an exact 
nonstatic solution of the 5D EYMGB theory in the pres- 
ence of a null fluid and by employing the Wu-Yang ansatz. 
We shall present a class of 5D nonstatic solutions describ- 
ing the exterior of radiating black holes with null fluid 
endowed with gauge charge, i.e., an exact Vaidya-like 
solution in 5D EYMGB theory. The Vaidya geometry 
permitting the incorporations of the effects of null fluid 
offers a more realistic background than static geometries, 
where all back reaction is ignored. It may be noted one of 
few nonstatic black hole solutions is Vaidya pLE] which is a 
solution of Einstein's equations with spherical symmetry 
for a null fluid (radially propagating radiation) source. It 
is possible to model the radiating star by matching them 
to exterior Vaidya spacetime (see [13, [3| for reviews on 
Vaidya solution and for it's higher-dimensional ver- 
sion). This Letter also examines the effect of the GB 
terms and YM gauge charge on the structure and loca- 
tion of the horizons for the radiating black holes. A black 
hole has three horizon like surface [ig, [i3| '■ time- like limit 
surface (TLS), apparent horizons (AH) and event hori- 
zons (EH). In general the three horizon does not coincide 
and they are sensitive to small perturbation. For a clas- 
sical Schwarzschild black hole (which does not radiate), 
the three surfaces EH, AH, and TLS are all identical. 
Upon switching on the Hawking evaporation this degen- 
eracy is partially lifted even if the spherical symmetry 
stays. We have then AH=TLS, but the EH is different 
from AH=TLS. In particular, the AH is located inside 
the EH, the portion of spacetime between the two sur- 
faces forming the so-called quantum ergosphere. If we 
break spherical symmetry preserving stationarity (e.g., 
Kerr black hole), then AH=EH but EH ^ TLS. Here 
the ergosphere is the space between the horizon EH=AH 
and the TLS, usually called the static limit [lB|- In both 
cases particles and light signals can escape from within 
the ergosphere and reach infinity. The characteristics of 
EH and AH associated with black holes in 5D EYMGB 
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are also discussed. 



II. VAIDYA-LIKE SOLUTION IN 5D EYMGB 
THEORY 

We consider SO{4) gauge theory with structure con- 



stant C^f^, the YM fields F^'^' and the YM potential 



(a) 



A'f^ . The gauge potentials A)^' and the Yang-Mills fields 
F^J^-* are related through the equation 



•^(/3)(7) 



F, 



(") 



ab 



We note that the internal indices {a, /3, 7, ...} do not differ 
whether in covariant or contravariant form. The action 
which describes EMYGB theory in 5D reads 0| : 



la 



M 



{R + uj'Lgb) 



^ pia)p{a)ab 
a=l 



(2) 

Here, g = dei{gab) is the determinant of the metric ten- 
sor, R is the Ricci Scalar and w' — uj/2 with uj the cou- 
pling constant of the GB terms. This type of action is 
derived in the low-energy limit of heterotic superstring 
theory [l^. In that case, uj is regarded as the inverse 
string tension and positive definite, and we consider only 
the case with > in this paper. Expressed in terms 
of Eddington advanced time coordinate (ingoing coordi- 
nate) w, with the metric ansatz of 5D spherically sym- 
metric spacetime [l^ [l^ [2^ : 



^A{v,rff{v,r) dv^ + 2A{v,r) dv dr + r^dill, (3) 



where dr2§ 



de^ 



+ sin 9 sin cfP'dilP'. Here 



A is an arbitrary function of v and r and {a;°} — 
{w, r, 6, (p, We wish to find the general solution 
of the Einstein equation for the matter field given by 
Eq. (|13l) for the metric ([3]), which contains two arbitrary 
functions. It is the field equation Gf ^ that leads to 
A{v,r) = g{v) [l^, [l^. This could be absorbed by writ- 
ing dv = g{v)dv. Hence, without loss of generality, the 
metric ([3]) takes the form , 



ds^ = -/(u, r)dv'^ + 2dvdr + r^dnj. 
We introduce the Wu-Yang ansatz in 5D d, 0| as 



(4) 



(5) 



2 < i < 4, 
1 < J < « - 1, 
!<(«)< 6, 

where the super indices a are chosen according to the 
values of i and j in order Q- It is easy to see that 
for the metric |3]) , the YM matter field equations admit 
solution a = Q [1, 0|- The Wu-Yang solution appears 



highly nonlinear because of mixing between spacetime 
indices and gauge group indices. However, it is linear 
as expressed in the nonlinear gauge fields because purely 
magnetic gauge charge is chosen along with position de- 
pendent gauge field transformation jj. The YM field 
2-form is defined by the expression 

= + A A^^'l (6) 

The integrability conditions 

dfM + iC(%\^)A(^)AF(^)=0, (7) 

as well as the YM equations 

d * F^'^^ + A = 0, (8) 

are all satisfied. Here d is exterior derivative, A stands 
for wedge product and * represents Hodge duality. All 
these are in the usual exterior differential forms notation. 
The GB Lagrangian is of the form 

LgB = i?' - mabR''^ + RabcdR^''"''. (9) 

The action ([2]) leads to the following set of field equations: 



where 



Gab — Gab + uj' Hab = Tab 



Gab — Rab — -^QabR, 



(10) 



(11) 



is the Einstein tensor and 



Hab — 2[RRab ~ 2RaaRb — '^R"'^ Raabp + R"^^ Ri 



-■^gabLcB, 



■baP'yl 

(12) 



is the Lanczos tensor. The stress-energy tensor is written 
as 

Tab^T^b + T^b: (13) 

where the gauge stress-energy tensor is 



Tab - 



Op(a)Ap(a) 
^-^a ^bX 



1 



^ \a ^ yab 



(14) 



The energy-momentum tensor of a null fluid is 

Tab^^{v,r)l3aPb, (15) 

where V'(^^: is the non-zero energy density and /3a is a 
null vector with 



/3a = 5°,/3a/3'^=0. 



(16) 
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Introducing 

Xi — r cos?/; sin sin 0, 

X2 — rsini/;sin0sin0, 

a;3 = rcos</)sin0, 

X4 = rcos0, 

and using ansatz ^ one obtains 
yl^^-* = Q sin'^ 61 (cos ?/; d0 — cos (/)sin-!/) sin0 c?-)/)) . 



^(3) = Q 



sm 



' (sin ip dcj) + cos cj) cos ip sin ^ c?-0) . 



A'-^^ = Q (sin 6* (cos cos (j)d(j) — sin i/; sin 4> dip) cos 6* 

+ cos ip sin c?^) , 

^(5) ^ Q (cos (/) (i6' - cos 6* sin sin 6* cf0) , 

= Q (cos 0(i6' - cos 6* sin sin 6* d(/)) . 



We then find that the non-zero components would read 
as: T; = ip{v,r), = T; = -3QV(2r4) and = 
= = Q'^/r'^. It may be recalled that energy- 
momentum tensor (EMT) of a Type II fluid has a double 
null eigenvector, whereas an EMT of a Type I fluid has 
only one time-like eigenvector [l^, [2lj . It may be noted 
that, the gauge held has only the angular components, 
Fg'.g^ with i ^ j, non-zero and they go as which in 
turn make go as r^*. 

The only non-trivial components of the EGB tensor 
(G^), in a unit system with cj' = w/2, take the form: 



g: = g: 



4:UJ , 



/'--(l-/) + —(l -/)/', 



(17) 



Gl = Gt = Gt = /"+-/'+^(l-/)+^ [/"(I - /) + f 

(18) 



a: = 1^ + ^/(1-/). 
2 r r'^ 



(19) 



Then, f{v,r) is obtained by solving only the (ITUl) . the 
equation G^ = is integrated to give the general solu- 
tion as 



4ujM{v) 8wQ2inr^4w2 
(20)" 

where M{v) is positive and an arbitrary function of v 
identified as mass of the matter. The gauge charge Q 
can be either positive or negative. The special case in 
which M{v) = and = 0, Eq. ^ leads to GB- 
Schwarzschild solution, of which the global structure is 
presented in The solution pO|) is a general spheri- 
cally symmetric solution of the 5D EYMGB theory with 
the metric ^ for the null fluid defined by the energy- 
momentum tensor p^ . Since YM go as r""* (the 



same as for Maxwell field in £> = 4), for 5D. That is 
why its contribution in f{v,r) will be the same for 5D 
as in 4D Reissner-Nordstrom (RN) black hole The 
nonradiating limit of this would be 5D-Yaskin black hole 
and not 5D analogue of Reissner-Nordstrom. 

There are two families of solutions which correspond to 
the sign in front of the square root in Eq. ([20)1 . We call 
the family which has the minus (plus) sign the minus- 
(plus-l-) branch solution. From Gl = , we obtain the 
energy density of the null fluid as 



-!/'(«, r) 



3 M{v) 
2 r3 



(21) 



for both branches, where the dot denotes the derivative 
with respect to v. We flrst turn our attention to the three 
limiting cases when the solution is known. These are (i) 
fil: Q = and a; 7^ then 5D-EGB black holes 
123, [ill. The solution of the Eq. (HH) is 



4ujM (v) 



(22) 



(ii) M(v) ^0, Q^Oandw = then the 5D-EYM black 
holes [20] ■ Now one has solution of the Eq. p7)) as 



f{v,r) = l- 



M{v) 2Q2lnr 



and (iii) in the general relativistic limit uj 
0, the minus-branch solution reduces to 



f{v,r) = 1 



M{v) 



(23) 

and Q2 ^ 
(24) 



which is the 51? Vaidya solution [l^ [T^ in Einstein the- 
ory. It may be noted that, in 513 Einstein theory, the 
density is still given by Eq. (|2ip. There is no such limit 
for the plus-branch solution. The family of solutions dis- 
, cussed here belongs to Type II fluid. However, In the 
static case M = constant and the matter field degener- 
ates to Type I fluid, we can generate static black hole 
solutions obtained in [1, 0] by proper choice of these con- 
stants. In the static limit, this metric can be obtained 
from the metric in the usual spherically symmetric form. 



with 



I{r) 



ds' 



r 

1+2^ 



-fir) df 



dr^ 
W) 



r'idns)^. 



(25) 



1±A 1 



4ujM 8ujQ^\nr Acu^ 

r.4 ^ r.4 



(26) 

if _j. this solution reduces to the solution which was 
independently discovered by Boulware and Deser [l[ and 
Wheeler {2:|. 

The Kretschmann scalar {K = Rabcd.R'^^'^'^ , Rabcd is 
the 5D Riemann tensor ) and Ricci scalar [R = RabR"^'', 
Rab is the 5D Ricci tensor) for the metric (UJ reduces to 



K = f" 



(27) 
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and 



where, 



i? -/" + -/' -4(1-/)' 
r 



(28) 



Radial {9 and — const.) null geodesies of the metric 
(|4]) must satisfy the null condition 



dr 



1±W1 + 



4wAf(u) 



8cjg2 1nr 4a;2 



(29) 

The invariants are regular everywhere except at the ori- 
gin r = 0, where they diverge. Hence, the spacetime 
has the scalar polynomial singularity [2l| at r = 0. The 
nature (a naked singularity or a black hole) of the singu- 
larity can be characterized by the existence of radial null 
geodesies emerging from the singularity. The singularity 
is at least locally naked if there exist such geodesies, and 
if no such geodesies exist, it is a black hole. The study 
of causal structure of the spacetime is beyond the scope 
of this paper and will be discussed elsewhere ^25!] . 
Energy conditions: The family of solutions discussed 
here, in general, belongs to Type II fluid defined in [2l|. 
When m — m{r), we have 'i/'=0, and the matter field 
degenerates to type I fluid [Ij, [ij]. In the rest frame 
associated with the observer, the energy density of the 
matter will be given by 



M = r;, p = -tI = -t: 



(30) 



and the principal pressures are Pi = (no sum con- 
vention) and due to isotropy P — Pi for all i. a) The 
weak energy conditions (WEC): The energy-momentum 
tensor obeys inequality Tabw'^w^ > for any time-like 
^Z- > 0, p > 0, P > 0,. We say that 



vector 
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I.e. 



strong energy condition (SEC), holds for Type II fluid if, 
WEC is true., i.e., both WEC and SEC, for a Type II 
fluid, are identical [T3| . 

h) The dominant energy conditions (DEC) : For any 
time-like vector Wa, T"-'''Wa'Wb > 0, and T°'''wa is non- 
space-hke vector , i.e., ip > 0, p > P > 0- Hence WEC 
and SEC are satisfied if M (v) > 0. In addition DEC also 
holds. 



III. HORIZONS OF 5D RADIATING BLACK 
HOLE 

The line element of the radiating black hole in 5D 
EYMGB theory has the form ([4]) with f{v,r) given by 
Eq. ||^D|) and the energy-momentum tensor p^ . The lu- 
minosity due to loss of mass is given by L = —dM/dv, 
where L < 1. Both are measured in the region where 
d/dv is time- like. In order to further discuss the physical 
nature of our solutions, we introduce their kinematical 
parameters. As first demonstrated by York [l^ and later 
by others the horizons may be obtained to 0{L) 

by noting that a null- vector decomposition of the metric 
^ is made of the form 



gab 



-PJb ~ la(3b + lab, 



(31) 



Pa 
lab 



1 



f{v,r)S: + 5l 



r'5X + r^sin\e)5t5t 
+r^sin\e)^in\c^)5t8t, 

lab /S'^O, 



(32) 

(33) 
(34) 



with f{v, r) given by Eq. ([20)) . The Raychaudhuri equa- 
tion of null-geodesic congruence is 



^ = «e - Rabi^i' - ^e^ 

dv 2 



ab 



(35) 



with expansion 0, twist $7, shear cr, and surface gravity 
K. The expansion of the null rays parameterized by v is 
given by 



(36) 



where the V is the covariant derivative and the surface 
gravity is 

K=-p''l''Vbla- (37) 

In the case of spherical symmetric, the vorticity and shear 
of la are zero. Substituting Eqs. (gHl) and into ((57| . 
we obtain surface gravity 



r 
+ 



AujM{v) 8a;Q2lnr 4cj2 



38) 



2M(v) q2 4Q^ Inr , 2lo 



y-L H ipi ;3 H -pr 



Then Eqs. dSO]), (|32]), ([361), and (09]) yield the expansion 
of null ray congruence: 



r 

4w 



4cjAf(w) 8wQ2inr 4^2 



The apparent horizon (AH) is the outermost marginally 
trapped surface for the outgoing photons. The AH can 
be either null or space-like, that is, it can 'move' causally 
or acausally [T6| . The apparent horizons are defined as 
surface such that ~ which implies that / = 0. From 
the Eq. (l39l) it is clear that AH is the solution of 



1-f- 
2w 



1 



1 



AloM{v) 8wg2 1nr 4a;2 
13 I4 ^ ,,4 



i.e., zeros of 



r"^ - M{v) + 2Q^\n{r) = 0. 



= 0. 

(40) 

(41) 



For Q and constant M, we have 5D Schwarzschild 
horizon r = ^/M. In general, Eq. (|41l) . which admits 
solutions 



riAH = exp 



1 Q2LambertW (0, x) + M{v) 
2 



(42) 



5 



To AH = exp 



1 g^LambertW (-1, x) + M{v) 
'2 



(43) 



Herex = -2SE(^^|(H)Zfi!). 

Here riAH and tqah are respectively inner and outer 
horizons and the LambertW function satisfies 

LanibertW(2:) exp [LanibertW(a;)] = x. 

The important feature of Eq. (|43l) is that it is w inde- 
pendent. This lead to fact that it is similar to pure 
EYM case. Thus the GB term does not cause the AHs 
of the 5D EYM black holes to be distorted. The TLS 
for a black hole, with a small luminosity, is locus where 
gvv — 0. Here one sees that 9 = 0, implies / = or 
gwii" = ^ah) = implies that r = tah is TLS and AH 
and TLS coincide in our non-rotational case. The pure 
charged case {M(v) = 0) is also important, then we have 
horizon without mass 



riAH = exp 



1 / 1 
-LambertW 0, --r 

2 \ 



To AH = exp 



-^LambertW ( — 1, 7^ 



(44) 



(45) 



For an outgoing null geodesic, r is given by Eq. ([291 
Differentiation of ([^^ w.r.t. v gives 



2w 



1 - \ 1 



4wM(w) 



8^(32 In r 4^2 



;46) 



L I 1M{y)r 



4Q Inrf 
— ;3 — 



2wr 



1 



3 



8wQ2 Inr 
3 



At the time-like surface r = tah, = and r > for 
L > 0. Hence photons escape from r = tah to the to 
reach arbitrarily large distances from the hole. 

However, In the general GB term does change the lo- 
cation of AH, e.g., in the limit Q — >■ , in the 5D-EGB 
case the AHs reads [l9ll 



TAH 



= {v) - 2uj 



Further, In the relativistic limit 



(47) 

0, Q ^ then 
TAH ^jM[v). 

The future event horizon (EH) is the boundary of the 
causal past of future null infinity, and it represents the 
locus of outgoing future-directed null geodesic rays that 
never manage to reach arbitrarily large distances from 
the hole. This definition requires knowledge of the en- 
tire future history of the hole. However, York [l6|, for 
a radiating black holes, argued that the question of the 
escape versus trapping of null rays is, physically, matter 
of qualitative degree and proposed a working definition 
of definition as follows: the EH are strictly null and are 
defined to order of 0(V) and Photons are in captivity by 



event horizon for times long compared to dynamical scale 
of the hole. It can be seen to be equivalent to the require- 
ment that for acceleration of null-geodesic congruence at 
the EH, 



0. 



EH 



(48) 



This criterion enables us to distinguish the AH and 
the EH to necessary accuracy. An outgoing radial null 
geodesic which is parameterized by v satisfies 



dr _\ 
~ 2 



1 



r 



/ 4uM(v) 8a;Q2lnr 4^2 

(49)" 



Then Eqs. §9^ and (|39l) can be used to put Eq. (gHl) in 
the form 

3 df] 



EH 



2r dv 
1 



EH 



3L 



EH , 1 



2r 

where the expansion is 
Qeh 



iujM(v) _ 8ujQ^ Inr EH , 4w£ 



EH 



m 



2rEH 



4wM(w) 8a;Q2lnr£H 4^2 



' EH 



' EH 



' EH 



(51) 



For the null vectors la in Eq. (1321) and the component of 
energy-momentum tensor yields 

3 df 



Rabl'^l 



a ih 



2r dv 



(52) 



The Raychaudhuri equation, with a ^ Q = 

dv 2 

Since EH are defined to 0{L), we neg lect 92, as 92 
0(L2). Eqs. dSO]), (HI]) and imply that 



(53) 



dB_ 

dv 



0. 



(54) 



EH 



Following [16|, ll9| , for low luminosity, the surface grav- 
ity K can be evaluated at AH and the expression for 
the EH can be obtained to 0{L). It can be shown 
that the expression for the event horizon is the same as 
that for the apparent horizon with M being replaced by 
M* [1^, where M* is effective mass defined as follows: 
M*{v) = M{v) - Ljn. From the Eq. (gZl), it is clear 
that, in general, the presence of the coupling constant, 
of the GB terms, lo produces a change in the location of 
horizons. Such a change could have a significant effect in 
the dynamical evolution of these horizons. For nonzero 
the structure of the horizons is non-trivial. However, the 
eq. (|43)) is independent of the of the GB coupling constant 
w, i.e., AH are exactly same as that in EYM without GB 
coupling constant w. Thus the GB term does not alter 
the horizons of the 5D EYM black holes. 
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IV. DISCUSSION AND CONCLUSION 

In this Letter we have obtain an exact black hole so- 
lution that describes a null fluid in the framework of 
5D EYMGB theory by employing 5D curved-space gen- 
eralization Wu-Yang ansatz. Thus we have an explicit 
nonstatic radiating black hole solution of 5D EYMGB 
theory. We have used the solution to discuss the conse- 
quence of GB term and YM charge on the structure and 
location of the horizons 5D radiating black hole. The 
AHs are obtained exactly and EHs are obtained to first 
order in luminosity by method developed by York [l6l |. 
We shown that a 5D radiating black hole in EYGB has 
three important horizon-like loci that full characterizes 
its structure, viz. AH, EH and TLS and we have rela- 
tionship of the three surfaces veh < t^ah — ^tls and 
the region between the AH and EH is defined as quan- 
tum ergosphere. The presence of the coupling constant 
of the Gauss-Bonnet terms uj produces a change in the 
location of these horizons [l^ . Such a change could have 
a significant effect in the dynamical evolution of these 
horizons. However, It turns out that the presence of the 
coupling constant of the GB terms oj > does not al- 
ter the location of the horizons from the analogous EYM 
case, i.e., horizons of the 5D EYM and 5D EYMGB are 
absolutel y sa me when obtained by procedure suggested 
by York [ig] to 0{L) by a null- vector decomposition of 
the metric. 



In 4D, the Vaidya-like solution of EYM yields the same 
results as one would expect for the charged null fluid in 
EM theory, i.e., the geometry is precisely of the charged- 
Vaidya form and the charge that determines the geom- 
etry is magnetic charge. This is because go over 
which is exactly the same as energy momentum of 
EM theory. However, this does not hold in 5D case be- 
cause components of energy-momentum tensor for EM 
and EYM theories are not same. Thus the Yaaskin's Q 
theorem does not hold in 5D case. The 5D solution dis- 
cussed here incorporates a logarithmic term unprecedent 
in 4D. 

The family of solutions discussed here belongs to Type 
II fluid. However, if M = constant and the matter field 
degenerates to Type I fluid, we can generate static black 
hole solutions obtained in [y, 0| by proper choice of these 
constants. In particular, our results in the limit w — ^ 
and Q reduce vis-a-vis to 5D relativistic case. 
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